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Abstract. Galactic cosmic rays are believed to be accelerated at supernova remnants via 
diffusive shock acceleration. Though this mechanism gives fairly robust predictions for the 
spectrum of particles accelerated at the shock, the spectrum of the cosmic rays which are 
eventually injected in the interstellar medium is more uncertain and depends on the details of 
the process of particle escape from the shock. Knowing the spectral shape of these escaping 
particles is of crucial importance in order to assess the validity of the supernova remnant 
paradigm for cosmic ray origin. Moreover, after escaping from a supernova remnant, cosmic 
rays interact with the surrounding ambient gas and produce gamma rays in the vicinity of 
the remnant itself. The detection of this radiation can be used as an indirect proof of the fact 
that the supernova remnant was indeed accelerating cosmic rays in the past. 
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Supernova remnants (SNRs) are believed 
to accelerate galactic cosmic rays (CRs) up to 
at least the energy of the knee (« 3 PeV) via 
diffusive shock acceleration (e.g. Hillas 2005). 
Though this is a very popular idea, a conclu- 
sive proof of its validity is still missing. The 
detection of a number of SNRs in GeV and 
TeV gamma rays (e.g. Caprioli 2011) is en- 
couraging, since gamma rays are expected to 
be produced in hadronic interactions between 
the accelerated CRs and the ambient gas swept 
up by the shock. However, competing leptonic 
processes may also explain these observations, 
and thus the debate on the origin of such emis- 
sion is still open. Notably, the emission from 
the SNR RX J17 13.7-3946 seems to be lep- 
tonic (Abdo et al. 201 1), the one from Tycho 
hadronic (Giordano et al. 2011), while for 
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most of the other gamma ray bright SNRs the 
situation is less clear. 

Gamma rays can also be produced by CRs 
after they escape from SNRs due to hadronic 
interactions in the ambient gas. Thus, the de- 
tection of gamma ray emission from the vicini- 
ties of SNRs might provide an indirect evi- 
dence for the fact that the nearby SNR is (or 
was) accelerating CRs (Aharonian & Atoyan 
1996; Gabici & Aharonian 2007). Massive 
molecular clouds have been detected in TeV 
gamma rays in the vicinity of the SNR W28 
(Aharonian et al. 2008). This might constitute 
a first example of such situation (Gabici et al. 
2010). Future observations with the Cherenkov 
Telescope Array will undoubtedly increase the 
number of detected SNR/cloud associations 
and clarify the issue (CTA consortium 2010). 

Unfortunately, the way in which CRs es- 
cape SNRs is little understood. Qualitatively, 
CRs with the highest energies are expected to 
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escape first, while CRs of lower and lower en- 
ergy are expected to escape progressively later 
in time, as the shock slows down (Ptuskin & 
Zirakashvili 2003), but the details of how this 
happens are still unclear. In this paper, the 
present status of the studies on CR escape from 
SNRs is reviewed. A review of the aspects re- 
lated to the gamma ray radiation produced by 
runaway CRs in the surrounding of SNRs can 
be found in Casanova, these proceedings. 

Qualitative estimates 

A supernova explosion drives an expanding 
shock in the interstellar medium. As long as the 
mass of the material ejected in the explosion is 
much smaller that the mass of the gas swept 
up by the shock, the shock speed u s y, remains 
roughly constant. This stage of the SNR evo- 
lution is called free-expansion phase. If parti- 
cles are accelerated at the shock via diffusive 
shock acceleration (e.g. Drury 1983), the ac- 
celeration time is: t acc « D/u 2 h , where D is the 
particles diffusion coefficient. Assuming Bohm 
diffusion (D oc E/B, where E is the particle en- 
ergy and B the magnetic field strength assumed 
here to be constant in time), as seems appro- 
priate at shocks, and equating the acceleration 
time to the age of the SNR, t age , one gets the 
maximum energy that a particle can reach at 
the shock: E max oc t age . Thus, in this phase, the 
maximum energy increases linearly with time. 

As a consistency check, note that CRs with 
the highest energy diffuse ahead of the shock 
one diffusion length Id as D(E max )/u s h before 
returning to it. It is straightforward to show that 
this length is of the same order of the SNR ra- 
dius R s h = u s htage- Thus, by adopting a Hillas- 
like criterium for CR confinement in the SNR 
(diffusion length < SNR size) one can conclude 
that all the accelerated CRs are effectively con- 
fined and particle escape does not occur 1 . 

Once the mass of the gas swept up by the 
shock becomes larger than the mass of the 

1 CR escape may indeed occur also during the 
free expansion phase (Ptuskin & Zirakashvili 2005) 
due to the actual, though slow, decrease of the shock 
velocity with time. However, the number of particles 
involved in the shock acceleration is quite small in 
this phase, so in the following we will neglect them. 



ejecta the SNR enters the adiabatic, or Sedov 
phase. For an expansion in an uniform medium 
this phase is characterized by the scalings: 
R sh oc f 2 / 5 and u s h °c r 3 ^ 5 . Let us call E* nax 
the maximum energy of the accelerated par- 
ticles at the time f of the transition between 
the free-expansion and the Sedov phases. Such 
particles are confined within the SNR because 
their diffusion length is of the same order of 
magnitude than the SNR radius ld(E^ ax , f) as 
R s h(t*)- Once the SNR enters the Sedov phase, 
the particle diffusion length increases as Id as 
D{E* max )lu sh oc t 3/5 , faster than the SNR radius, 
which scales like f 2 ^ 5 . Therefore, particles with 
energy E* max , which were confined during the 
free expansion phase, now violate the Hillas 
criterium and can escape the SNR. This implies 
that the decrease of the shock velocity reduces 
the capability of the SNR to confine particles. 

Qualitatively, to see when CRs with differ- 
ent energies leave the SNR, one can equate the 
particle diffusion length with the SNR radius, 
to get: E max oc R sh u sh B oc r 1 ^ 5 . This very rough 
estimate has the sole purpose to provide us 
with a qualitative description of the the process 
of CR escape from SNRs. What can be under- 
stood from that is that CRs are released mainly 
during the Sedov phase, and that they are re- 
leased gradually, the ones having the highest 
energy first, and the ones with lower and lower 
energies at later and later times. The decrease 
with time of the maximum energy of confined 
particles is quite mild, but this has to be as- 
cribed uniquely to the fact that a diffusion co- 
efficient constant in time has been assumed in 
the derivation. If the magnetic field at the shock 
is amplified by CR streaming instability, the 
value of the diffusion coefficient is expected 
to increase as the shock decelerates (Ptuskin 
& Zirakashvili 2003). As a result, a faster de- 
crease of E max with time may be expected. 

The approach presented here is (deliber- 
ately) qualitative and suffers from many lim- 
itations. Nevertheless, it gives a fairly reason- 
able description of the main features and trends 
which are also derived from detailed models, 
which will be reviewed in the next Section. 

Another important issue related to the pro- 
cess of particle escape from SNRs concerns the 
shape of the CR spectrum which is eventually 
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released in the interstellar medium. It is well 
known from the measurements of the abun- 
dances of secondary elements such as Li, Be, 
and B in the cosmic radiation that the CR con- 
finement time in the Galaxy decreases with en- 
ergy as t esc oc E~ s , with s « 0.3. ..0.6. Thus, to 
reproduce the shape of the observed CR spec- 
trum oc E~ 2J ' , CR sources have to inject a spec- 
trum with slope close to, but definitely steeper 
than 2: N inj oc E~ 2J+S * £- 2 -i- 2 - 4 . How does 
this compare with the spectrum that SNRs are 
expected to inject in the interstellar medium? 

To estimate such a spectrum one can pro- 
ceed as follows (see Ptuskin & Zirakashvili 
2005, for a detailed treatment): let f esc (t) be 
the fraction of the kinetic energy flux across 
the SNR shock that is carried away by escap- 
ing particles. For a SNR expanding in a homo- 
geneous medium and in the Sedov phase, the 
energy carried away per unit time by escap- 
ing CRs reads: dS/dt = f esc (pu 3 J2)(4nR 2 sh ) oc 
f esc t~ l ■ Let us now assume that particles of en- 
ergy E escape at a time t oc E~ 5 . This allow 
us to write the spectrum of escaping CRs as: 
N esc = (d&/dt)/E x (dt/dE) oc f esc E~ 2 . Thus, 
if f esc is constant in time, the spectrum of the 
particles released by a SNR in the interstel- 
lar medium during the whole Sedov phase is 
oc E~ 2 , regardless of the value of the parameter 
6 which describes the way in which CRs of dif- 
ferent energies are released. Remarkably, such 
result has beed obtained without specifying the 
spectral shape of the particles at the shock. 

This last point deserves more attention 
(what follows is a simplified summary of the 
results derived by Ohira et al. 2010). Let us 
now assume, in a slightly different fashion, that 
a fraction rj (again, for simplicity, constant in 
time) of the kinetic energy flowing across the 
shock per unit time is converted into accel- 
erated particles. These particles have a spec- 
trum Q C R oc E~ a [#/time/energy] and after 
acceleration are advected downstream of the 
shock. The source term for escaping particles 
can be written as: Q esc oc E Qcr(E, t) S(E - 
Emax(t)), where we have assumed that only 
recent acceleration is relevant (Ptuskin & 
Zirakashvili 2005). In an expanding SNR par- 
ticles accelerated at earlier times are advected 
downstream and suffer adiabatic energy losses, 



and thus their contribution can, as a first ap- 
proximation, be neglected. The energy per unit 
time in form of CRs that escape the shock is: 
Fesc = J Qesc E AE oc Q CR (E ma ^)E 2 max . 

First, consider the situation in which the 
slope of the spectrum of the accelerated CRs 
a is smaller than 2. In this case, the particles 
with the highest energy carry most of the to- 
tal energy in form of CRs, and we can write: 
r]L k * QcR(E max )E 2 max oc Fesc- This means 
that a constant fraction of the shock kinetic 
energy flux is converted into escaping parti- 
cle, which, as discussed above, implies that the 
spectrum of the particles released by the SNR 
is Nesc x E~ 2 . Conversely, if the spectrum of 
the CRs accelerated at the shock is softer than 
2, particles with energy Eq ~ mc 2 carry most 
of the energy, so we have: tjL^ « Qcr(Eo)E 2 ) oc 
QcR(E max )E 2 lax (E /E max ) 2 - a . The last term in 
brackets tells us that f esc oc E 2 ~° x and therefore 
in this situation the spectrum of escaping CRs 
is softer than a = 2 and equal to N esc oc E~ a . 
Thus, in order to be consistent with CR data, 
not only the spectrum injected in the interstel- 
lar medium, but also the spectrum of CRs at 
the SNR shock has to be steeper than E~ 2 . This 
will be discussed in the next Section. 

Overview of models 

The acceleration of particles at shock waves 
is believed to be a very efficient mechanism, 
i.e. a large fraction of the shock kinetic energy 
flux can be converted into accelerated parti- 
cles. As a consequence, the pressure term due 
to the presence of accelerated particles can- 
not be neglected in the flow equations (as it 
is done in the test-particle limit), the shock 
structure is then modified by the presence of 
CRs, and the whole process of acceleration be- 
comes non-linear (for a review see Malkov & 
Drury 2001). Extensive efforts have been de- 
voted to the solution of the steady-state equa- 
tions describing non-linear shock acceleration, 
and an approximate analytic solution has been 
found by Malkov (1997). Accelerated par- 
ticles exhibit a concave spectrum, which is 
harder than E~ 2 (up to E~ l ,5 ) at the high- 
est energies. In Malkov's aproach, as well as 
in several other analytic or semi-analytic ap- 
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proaches, the maximum momentum of acceler- 
ated particles, p max , is a free parameter, and the 
particle distribution function vanishes abruptly 
above that momentum. Since a flux of particles 
upward in momentum exists at the shock, this 
is equivalent to assume that particles acceler- 
ated over p max escape the system (see Eichler 
1984, for an early discussion of this issue). 
Within this approximation, it is possible to esti- 
mate the amount of energy subtracted from the 
system by the runaway CRs. Due to the hard- 
ness of the spectrum, this can amount to a sig- 
nificant fraction of the total bulk energy flow- 
ing across the shock (e.g. Berezhko & Ellison 
1999; Blasi et al. 2005; Amato et al. 2008). 

Another possibility to model particle es- 
cape is to assume a spatial escape bound- 
ary, instead of the one in energy space de- 
scribed above. This is more justified, since 
the CRs themselves are believed to amplify 
the magnetic field at the shock and generate 
the turbulence they need to be scattered off 
via resonant (Lagage & Cesarsky 1983) or 
non-resonant (Bell 2004) streaming instabil- 
ity. Depending on the CR flux at a given po- 
sition, the level of the magnetic turbulence 
is expected to decrease with the distance up- 
stream of the shock. This suggests the exis- 
tence of a spatial boundary beyond which CRs 
cannot be effectively scattered and freely es- 
cape. Reville et al. (2009) studied the acceler- 
ation of CRs at a shock in presence of a spa- 
tial escape boundary and assumed the turbu- 
lence to be generated at such boundary by the 
non-resonant streaming instability. This insta- 
bility injects short- wavelength turbulence, not 
very efficient in scattering high energy parti- 
cles with a large gyroradius, but in its nonlin- 
ear evolution is characterized by the appear- 
ance of large scale structures like cavities (Bell 
2005) which might reduce the mean free path 
of high energy particles. Thus, Reville et al. 
(2009) made the crude, but not unreasonable, 
approximation that the particle diffusion coef- 
ficient is of the Bohm type in a background 
magnetic field amplified to the value at which 
the non-resonant instability is believed to satu- 
rate. Within this approximation, the spatial es- 
cape boundary is expected to be roughly lo- 
cated at the position at which the turbulence 



growth rate equals Uadv/lp, where u ac i v is the 
velocity at which waves are advected towards 
the shock and l p is the distance upstream of the 
shock. The non-resonant mode of streaming in- 
stability is likely to dominate over the resonant 
one at large distances from the shock (Pelletier 
et al. 2006), and thus its properties should fix 
the position of the boundary. 

Though a fully self consistent treatment of 
the problem is still not available, introducing 
a spatial escape boundary presents several ad- 
vantages. First, it allows to determine the shape 
of the cutoff in the CR spectrum and the spec- 
trum of the escaping particles, which are no 
longer an abrupt step-function feature and a 
delta function in energy, respectively, as in the 
case of an energy escape boundary. This is very 
important in order to investigate the radiative 
signatures related to CR acceleration at shocks, 
whereas it is of minor relevance for the esti- 
mate of the radiation from runaway particles, 
whose spectrum is found to be quite narrow 
and peaked around p max (Reville et al. 2009; 
Caprioli et al. 2010b). Moreover, it is possible 
to check, a posteriori, the consistency between 
the assumed position of the escape bound- 
ary and the value of the amplified magnetic 
field/diffusion coefficient, that depends on the 
spatial CR intensity (Reville et al. 2009). 

Free escape boundaries have been used in 
Monte Carlo approaches also, both in energy 
(e.g. Ellison & Eichler 1984), and space (e.g. 
Vladimirov et al. 2006), as well as in time- 
dependent numerical simulations of shocks 
(e.g. Kang & Jones 1995), with the value of 
the maximum energy or the spatial position of 
the boundary as a free parameter of the models. 
A critical discussion on the use of free escape 
boundaries can be found in Drury (201 1). 

Zirakashvili & Ptuskin (2008b) also stud- 
ied the acceleration of CRs at a plane shock 
in presence of the non-resonant streaming in- 
stability and computed the CR diffusion coef- 
ficient in the resulting short-wavelength turbu- 
lence. A free escape boundary was imposed at 
a distance L upstream of the shock, with L of 
the order of the radius R s of the SNR to be 
simulated. One of the results from this simu- 
lation is that the width of the CR distribution 
upstream of the shock is small, of the order of 
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< 0.1 L, which implies that the CR acceleration 
at SNRs can be described with satisfactory ac- 
curacy as a one-dimensional phenomenon even 
in three-dimensional systems such as SNRs. 
Moreover, the maximum energy of accelerated 
CRs was found to be lower than the one that 
would be obtained by assuming Bohm diffu- 
sion in the amplified magnetic field, since the 
scattering in the short-wavelength turbulence 
is not so efficient. In particular, only type Ib/c 
and lib supernovae were found able to accel- 
erate CRs up to PeV energies, while for type 
II-P supernovae, the ones with the highest ex- 
plosion rate, a maximum energy of only a few 
hundreds of TeV was obtained. However, in or- 
der to fully understand how particles with the 
highest energy are confined and accelerated at 
a shock, and thus how SNR can accelerate par- 
ticles up to the PeV range and above, a study of 
the instabilities that involve long-wavelengths 
(i.e. larger than the gyroradius of the generat- 
ing CRs) is needed. Promising results in this 
direction have been recently published (Bykov 
et al. 201 1; Shure & Bell 201 1). 

To conclude this section we discuss, in the 
light of the detailed models described above, 
which are the expectations for the total CR 
spectrum injected by a generic SNR in the 
interstellar medium. As discussed in the pre- 
vious section, in order to be consistent with 
CR data, CR sources have to inject spectra 
steeper than E~ 2 , with slopes in the range « 
2.1 ...2.4. However, also the observed high level 
of isotropy of CRs has to be taken into ac- 
count. For example, an injection spectrum oc 
E~ 2A corresponds to a CR escape time from 
the Galaxy that increase quickly with energy 
t esc oc E - 6 . In this case, the escape time of 
high energy particles would be so short to re- 
sult in a strong anisotropy of CR at the Earth, in 
contrast with observations (e.g. Hillas 2005). 
Though the discussion on this is still ongoing 
and requires further investigation, let us choose 
as a reference value for what follows a slope of 
the injection spectrum equal to « 2.4, which 
should minimize the problem of anisotropy. 

Such a steep spectrum for the runaway 
CRs corresponds, roughly, to an equally steep 
spectrum for the CRs accelerated at the SNR 
shock (see previous Section). This seems to be 



in disagreement with the predictions of non- 
linear shock acceleration theories, according to 
which spectra are, at high energies, harder than 
E~ 2 . A way to solve this problem is the fol- 
lowing (Zirakashvili & Ptuskin 2008a): the 
spectrum of the accelerated particles is deter- 
mined by the compression factor of the veloc- 
ities of the scattering centres, which can be 
different from the shock compression factor r 
(ratio between up- and down-stream gas ve- 
locities), which is often used to determine the 
CR spectrum. If the scattering centres move 
against the fluid upstream and with it down- 
stream (as might be the case for Alfven waves), 
the compression factor that CRs feel at the 
shock may be reduced with respect to r and the 
resulting CR spectrum significantly steepened. 

By including this effect, named Alfven 
drift, in a numerical code describing non-linear 
CR acceleration at a spherical SNR shock, 
Ptuskin et al. (2010) obtained a good fit to the 
entire CR spectrum up to an energy of ~ 3 x 
10 18 eV, well above the knee. They considered 
the contributions from different types of super- 
novae (la, IIP, Ib/c, lib) and parametrized the 
value of the magnetic field at the shock in order 
to mimic the magnetic field amplification, as it 
is observed from X-ray synchrotron filaments 
(Volk et al. 2005). Bohm diffusion in the am- 
plified field was adopted. Less optimistic con- 
clusions have been drawn by Caprioli et al. 
(2010a), which employed a (one-dimensional 
and quasi-steady-state) model for non-linear 
particle acceleration to compute the spectrum 
of CRs escaping from individual SNRs, though 
they considered the magnetic field amplifi- 
cation due to resonant streaming instability 
(Lagage & Cesarsky 1983). They obtained 
power law spectra roughly of the order of E~ 2 , 
too hard to explain the observed CR spectrum. 
When including the effect of Alfven drift, spec- 
tra become steeper but do not reach the energy 
of the knee. Additional (non-resonant?) field 
amplification might mitigate this problem. 

Gamma rays from runaway CRs 

After escaping the SNR, CRs interact with the 
ambient gas and produce gamma rays. The 
production of such radiation is enhanced if 
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dense gas (e.g. a molecular cloud) is present in 
the vicinity of the SNR. The CR background 
in the Galaxy has a steep spectrum as 6 x 
lCr 3 (£7Te\ r r ' 7 eV crrT 3 , and thus an excess 
in the CR intensity would appear more easily 
at high energies. For this reason we focus here 
on the energy domain probed by Cherenkov in- 
struments (> 100 GeV). 

To maintain the CR intensity at the ob- 
served level, each SNR has to provide « 
10 50 £ 50 erg in form of CRs («10% of the 
explosion energy). For a CR spectrum with 
slope 2.4 this corresponds to: E 2 Ncr(E) as 
2xl0 60 £ 5 o(£/TeV)- - 4 eV. If CR diffusion pro- 
ceeds isotropically, particles with a given en- 
ergy at a given time after escape are distributed 
roughly uniformly in a spherical region of ra- 
dius R» as V6 D t surrounding the SNR. The 
energy density of runaway CRs exceeds the 
one in the background if /?» < 130 E l ^pc. The 
average diffusion coefficient of TeV CRs in the 
Galaxy is as 10 29 cm 2 /s and thus the excess of 
CRs lasts for as R 2 J6D as 10 4 yr. Therefore, the 
search for gamma ray emission from molecular 
clouds located close to SNRs has to be focused 
on regions of size 100 - 200 pc around SNRs 
not much older than as 10 4 yrs. 

The detection of such radiation can serve 
as an indirect proof for the fact that SNRs are 
indeed the sources of CRs. Moreover, the fig- 
ures given above depend on the local (close 
to the SNR) value of the diffusion coefficient, 
which might differ significantly from the av- 
erage galactic one (e.g. due to enhanced CR 
streaming instability in the vicinity of a CR 
source). Thus, the characteristics of the radi- 
ation may also be used to constrain the prop- 
erties of CR diffusion (see Gabici et al. 2009 
and Casanova, these proceedings for details). 
To date, very few SNR/cloud associations have 
been detected in TeV gamma rays (e.g. the 
SNR W28), but future observations with the 
Cherenkov Telescope Array promise to in- 
crease the number of detected objects and pro- 
vide a test for this scenario. 
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